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ABSTRACT 


The radial diffusion equation with synchrotron losses Is 

solved by the Laplace transform method for near*equator ial ly 

mirroring relativistic electrons. The evolution of a power 

law distribution function is found and the characteristics of 

synchrotron burn -off are stated in terms of explicit parameters 

for an arfa i trafy •diffu^ coefficient of the form Dj^^ = D^L®. 

The peaking of the 10.4 cm volume emissivity from Jupiter at 

provides an estimate 0 q=“9x 1O*’* sec ' in the radiation 

belts; this value is suggested as the appropriate modification 

for an equatorial field strength of 4,2 Gauss of Birmingham 

et al,*s (1974) result,- Non-synGhrotron losses are included 

phenomenoi ogical I y and from the phase space densities reported 

by Mcllwain and Fill ius (19755 particle lifetime is esti- 

8 2 - 0 . 

mated as t^6x 10 L sec. Asymptotic forms for the distri- 
bution in the strong synGhrotron loss regime are provided. 
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INTRODUCT ION 


Prior to the Pioneer 10 and 11 missions to Jupiter a con- 
siderable effort was made to estimate the fluxes of energetic 
particles in the radiation belts. The basic starting point 
was the two-dimensional interferometer maps of non- thermal syn- 
chrotron radiation provided by Oerge (1966) and Branson (1968) 
and deconvolved by Beard and Luthey (1973) into equatorial emis- 
slvities. The radial diffusion equation with synchrotron losses 
was then solved (Birmingham et al., 197^; Coroniti, 197^j Stansberry 
and White, 1974) and parameters varied until a best fit with 
the equatorial emissivity was found. This procedure gave a 
determination of the radial diffusion coefficient which was 
in reasonable agreement with theoretical descriptions of the 
diffusion process based on atmospheric neutral winds driving 
a fluctuating dynamo electric field (Brice and McDonough, 1973; 
Coroniti, 1974). • 

The data provided by Pioneers lO and 11 have been consis- 
tent with the model of an intense core of electrons with pitch 
angles near 9CP confined to equatorial magnetic latitudes and 
undergoing inward radial diffusion. The pancaked nature of 
the pitch angle distribution was demonstrated by Van Allen et 
al, ( 1974 ), who modeled the distribution as facsin^Q and found 
that for the 21, 31 MeV detectors M 3.5, 4 gave good closure 
for the Inbound/outbound Pioneer 10 data between 4<L<12. 

The signature of inward radial diffusion is manifest in all 
flux diagrams, J(L), especially in the vicinity of the Jovian 
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sateHites; if measured phase space densities at constant first 
invariant decrease towards the planets, as in Fig. 1, an inward 
diffusive current must result when the particles' third invariant 
is violated. 

The purpose of this paper is to re-examine the radial dif- 
fusion problem in light of the measurements of Pioneers 10 and 
11. Previous analyses have used'models for the electron dis- 
tribution which differed significantly in form from those measured 
and all conclusions heretofore have relied on numerical computa- 
tions, with underlying physical relationships not manifest (viz., 
how are Birmingham et al.'s (1974) results modified by an improved 
measurement of the magnetic field strength?). The mounting 
evidence that radial diffusion is a fundamental magnetospheric 
transport process warrants an analytic solution with conclusions 
of general applicab? 1 1 ty. * 

We solve the steady-state radial diffusion equation with 
synchrotron losses for near-equatorially mirroring relativistic 
electrons, driven by a diffusion coefficient modeled as 

o 

by the Laplace transform method previously utilized by Goroniti 
(1974) . This procedure brings the analysis in close analogy 
with the problem of heat conduction in solids {Carslaw and Jaeger, 
1959) . Our solution is similar to that of Birmingham et al. 

(1974) , who used an eigenftmction expansion technique with a 
delta f'unction source of particles in the outer belt. We, how- 
ever, concentrate on the evolution of a distribution function 
of particles {a power law in momenta) whose properties are now 
fairly well known. The solution is shown to vary as the loss-free 
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solution up to a point where synchrotron losses become domi- 
nant. A precise statement of this synchrotron ''burn-off regime 
is found, where significant reduction of the phase space den- 
sity occurs due to rapid synchrotron degradation on a time scale 
comparable to that of radial diffusion. The Laplace transform 
technique has the advantage of providing an integral represents- 
tion for the solution from which asymptotic forms for the dis- 
tribution in the "bujm-off" regime can be obtained. For application 
to Jupiter, the peak of the 10.4 cm voliime emissivity at L= 1.8 
(Birmingham et al., 1974) provides a lower bound to the magni- 
tude of the diffusion coefficient (i.e. , diffusion is strong 
enough to transport electrons into at least L = 1.8). The decay 
of the volume emissivity at smaller L values is presumed to 
be due to synchrotron degradation of the distribution which 
provides an upper boxand to the diffusion coefficient. 

Since Mcllwain and Fillius (1975) have found significant 
decreases in the phase space density inside of lo which are 
non- synchrotron associated, we have included a phenomenological 
loss termi and and resolved the radial diffusion equation excluding 
synchrotron losses. It is argued that for a particle lifetime 
which varies as t = solution is a power lay liv the 

loss region, consistent with observation. From the slope of 
measured phase space densities, the magnitude of particle life- 
time, T^, is estimated. 

o . . - . ■ 

2*"Ct 

Finally, if the particle lifetime persists as.x = 
into the synchrotron emission region (LAi 2) , a complete solution 
which includes both phenomenological and synchrotron losses 


it 


II SOLUTION OF THE RADIAL DIFFUSION EQUATION WITH SYNCHROTRON 
LOSSES 


The basic equation governing the steady-state radial dif- 
fusion with synchrotron losses is (Schulz and Lanzerotti, 1974) 


P» t) 


_6_|Il 0 - 

6Lr? ‘^LL 6LJ 


^[<M>f]-^[<J)f] = 0 
- ( 1 ) 

f(x/ p, t) is the phase space distribution function averaged 
over phase angles conjugate to M, J, $ and n(x, t) = /dpf is 
the local density. M = first adiabatic invariant/ 

0 is the equatorial pitch-angle/ J = i5p,,ds is the second invariant, 
and L - r/Rj is the magnetic shell parauneter for a dipole field 

The quantities <M> and <J> are the bounce-averaged 
rates of change of M and J due to synchrotron energy loss. 

Throughout this paper we make the assumptions that most of the 
electrons mirror near the equator (9 ^it/ 2) and have ultrarelativistic 

■i 


B(L) = B„L 
O 


■J 

-3 


energies y “ 


(i -vYc^y 


»1. For this case, the synchrotron 


degradation of the second invariant is small since <J)<xcos 0. 

• 2 

and 6/bJ<J)«cos 9 (Coroniti, 1974); hence the last term In 
(1) will be dropped compared to <M>, which Is finite at 9 = tt/2. 

mc2 2 

We approximate M as — bXlT* 

The rate at which M changes due to synchrotron emission 
is, to lov/<?.st order, 


<M> 


4 e^B^(L) 

m c 


yM 


m c 


1 


15/2 




3/2 


( 2 ) 



If 0|^j^ ss equation (1) Is of the form 


( 3 ) 


a) Initial Value Problem 

The consequences of the synchrotron term can be understood 

if we drop the diffusion term in equation (3) and solve an initial 

value problem with f(t = 0) = f^(M) = 0 (11 - M) . 0 (x) 

o o 

is the Heaviside step function which introduces an arbitrary 
cutoff, M^/ at large values of the first invariant. N is the 
spectral index of the distribution and A is a normalization. 

Two cases are of interest for comparison with the diffusion 
problem developed in the next section: if N = 2 the solution 

is 


. r 

f.(M, N = 2, t) = A l- ■ 

IT 


0 ( M^ - M ) 


(4) 


where 


M. 


“o = 


(1 + 




If t » tg “ 


2L 


15/2 


2L 


15/2 


s.t)2 


(5) 


o°s 


(the synchrotron energy half-life of the 


cutoff) a significant burn-off of particles occurs at even 
though M < <- . The second ■ terra in equation ( 4 ) is a synchro- 
tron correction to the N *= 2.spectr\ara. For' N = 1 the solution 
is ' 
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( 6 ) 

with the same as (5) . The synchrotron correction is absent 
since is a "stationary" solution. 


f(M, N = 1, t) = 



.372-® <“o - “> 


b ) Dfffusive Boundary Problem 

The diffusive boundary value problem has much the same 
effect. If the diffusion process cannot transport particles 
in faster than their synchrotron degradation, burn-off should 
occur. If we estimate characteristic times from equation (3) 

• 'l'5'2 ' I 

■^SYNC” Oj/M ' ■'■R.'D. ~p |_a-2 

then 

■^SYNC " p^a+l iTz 

.... 

where 





( 8 ) 


roughly marks the transition from the synchrotron burn-off regime 
to the pure radial diffusion regime (weak synchrotron losses). 


In the steady-state equation (3) becomes 





19/2 J_ .a-2 3f 
dL^- 3L 




/ I ‘ ® 

Lett'ng 6 = y « 7^, then If G(s)s [ e”®^G(y)dy, 

0 


l2 Jl 4+ (a - 


6L 


no) , . 


where we require 


Urn e"®^G(L, y) 
y -• 0 


1 im, e M) 

M -* flo' 


0 ni) 


The solution of equation (10) subject to the boundary conditions 
G^L ® Lj) “ 0 at the Surface, Lj « 1, and G^L = 


/ 1 1 r»“+® 

m^/2 2ttI 


I (xr) I (xb) - I (xb ) I \xr )^ 

. e®YG (y). V. ■■ :v. V V, . .“v, ,l ( 12 ) 

* o' ^ I (xa ) I' (xb I - I (xb ) I (xa ) 

^ V “V V- ■ " V 


where tha contour lies to. the right of all poles of the Inte- 
grand, l_^(xr) Is the modified Bessel function and 


la + 11/21'/ D 


r = i,'2(a + M/2) ^ ,^-i(a + 11/2) ^ ^ l_-i(a + n/2) 

.1 o 

(13) 


i-£gI:} g? vo- -\ 


At L- we assume that the distribution can be represented by 


a power law f 


with 


o ’ JpTN+zT* H2) then becomes 


f(t. M) = 


dte‘ 


rr (xr) l_ (Xb)-I (xb) I (xr)1 


2/Mt 


3 ■ a 


(a + ll/2)^0 


a + 11/2 


(15) 


If V is an Integer the bracket is replaced by Its limiting value 

V/hile l+^(Xr) has a branch point at the origin, the com- 
bination In the denominator of the bracket Is a single-valued 
even function of Xj as a function of t, all of the poles of 
the bracket lie on the negative real axis (Gray and Mathews, 
1922 ). If N is an integer the contour can be closed by a semi- 
circle to the left and the integral is evaluated by the theory 
of residues in Appendix A, 

The case of N * 1 is relevant for synchrotron-emitting 
particles at Jupiter since the time-averaged power spectrum 
of non-thermal emissions is flat between 200 MHz < f < 3000 MHz 
(Carr and Gulkls, 1969); if the E - M spectrum var ies as f”*, 

6 is related to the spectral index fay 5 .= ^(N- 1) (see equation 
(36)). Figure 1 is a plot of phase-space densities reported 
by Mcllwain and Fillius (1975) from Pioneer 10 measurements 
of energetic electrons. For L> 10 the particles have a spectrum 
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which has N'\'3.5. Inside Europa and lo, the spectrxim of low 
energy particles (y - 2 - 16) hardens considerably anv< at L = 3 
has settled into a profile of N 1 ; the corresponding range 
of M is 6.4- 360 MeV/G. 

From Appendix A have for N = 2 


f(L, M, N = 2) 






-- Lun 


n=l 


'••Tv n®) 

vn^ViH / v\ v,n / 


(16) 

and for N = 1, neglecting edge effects at 


/ L ( 3"Ct ) ^ 

f(L. M. N = 1) = a(^) ^ 


A ^ y 

rt 


2 

.Lun 


n=l 


ttJ I 

\j 

1 Y 

\Ly) . n_ ; 

I'j fy n^') 


H j>) 

’ V f n ) 

1 - J 

V 



• WVv.n'’) -'-»)]} 

(17) 


-2 

where \ » 


2/^ 


^ and A , A- are functions of L defined 
Z n O 1 


(a + i1/2)"D 
in Appendix A; are defined such tha 
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'’v(%,n'*) \(Yv,n®) • solutions 

can be compared v/ith Birmingham's (1974) equation (11) for the 

Invar! ant- space distribution function. The similarity with 

_2 

the initial value problem is clear. When X -»■ 0 only the 

term survives to reproduce the loss-free solution to equation 

(9). The term jX is the synchrotron correction to the N = 2 

2 

spectrum. The series of Bessel functions contributes when X » 

and strong synchrotron degradation occurs. In fact, if -^=0, 

X'^ 

the series in equation (17) is a representation of -A^ of the 

type described by Erddlyi et al. (1953) . Weak synchrotron losses 

can be expected when the argument of the exponential is not 
—2 2 

small and ^ 3 ^ • A ;few of the values of ^ have been tabu- 

lated by Jahnke and Emde (1945) . In particular, if v = 1/2 
and b » 1, Yl = -- I -, g nr for a = « 1. For other 

values of v, , is close to ir. We conclude that provided 



2/m 12 

;(o + 11/2)2d'^^ 


(18) 


weak synchrotron losses prevail throughout the diffusion to 
the surface. 

I 

The more interesting case involves particles which experience 
weak losses at but undergo strong degradation at L. Inspection 
of the bracket of equation (14) reveals that the criterion for 
weak (strong) losses at L depends on whether ,{Xr) is in the 
•small argument (asymptotic form) for most of that part of the 
integration up to where the exponential cuts off at t = -1. 

In this case an asymptotic expression for the distribution function 


in the burn-off regime can be obtained by expanding the rele- 
vant Bessel functions in equation. (14) and then inverting the 

transform. This procedure is carried out in Appendix B. Thus, 

— 2“2 1 2 

the appropriate modification of (18) is X r and provided 

that 

(a + 11/2^0 

weak synchrotron losses will prevail from L to L^.. This result 

is similar to the estimate given by equation (7) but differs 

-2 

significantly by the factor of (a + 11/2) . We have arbitrarily 

_2 2 12 

chosen the transition criterion at X r inspection 

of the argument of the exponential in equations (B.6) or (B.8) 

—2 2 

suggests an alternative X r < 4 which modifies the result only 
slightly. Although the solution {A. 14) and in particular equa- 
tions (16) and (17) are valid for integer values of N, the burn- 
off criterion (19) and the asymptotic forms are valid generally. 

c) The Diffusion Coefficient 

The synchrotron emissions (which peak at for 10.4 cm 

radiation) span the frequencies Af = 200 - 3000+ MHz. I f the 
electrons emit predominantly at the characteristic frequency 
^c ” ° ^ for a surface magnetic field strength 

of 4.2 Gauss, Ay = 8- 32 at L = 1.8, corresponding to 
am - 22 r 352 MeV/G. Equation (l9 ) can be expressed as 
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( 3 . 1 X 10~^)B^^VMlMeV7G~y 
(a + 11/2)^ L® ^ 

1 ,7 X \or^^^y^Jf\{mz] ^ 

' (a + ll/2)2 L^+5/2 ) 


( 20 ) 


-The Pioneer 10 data (Mcllwain and Fillius, 1975) indicate that 
for AM = 6,4 - 360 MeV/G the spectral index N»1 at L = 3. The 
flat synchrotron spectrum is consistent wi th N = 1 throughout 
the strong synchrotron radiation region, which extends to at 
least L = 1.8, The flat sytichrotron spestrum suggests that 
the radiating electrons diffuse to L « 1.8 without a signifi- 
cant synchrotron degradation of the N = 1 spectrum, we use (20) 
with the upper value of M =352 MeV/G = 3000 HHz^ to arrive 
at an estimate cf the m'/niroum diffusion coefficient 0^ which 
is consistent with the preservation of the N » 1 spectrum. 

We find that the minimum value 0 . x 10^® = 0,88,- 2,0, 4,6, 

11 sec"^ for a = 4, 3 , 2, 1, Note that for the 3000 MHz emitting 
particles to diffuse to the surface with weak losses requires 

D^x 10^*^ a40, 50 , 64, 86. ~ " ' 

o • 


Ill 


NON SYNCHROTRON LOSSES 


a ) General Scaling Laws 


For weak synchrotron losses we need retain only the 
term in the distribution -Function, and since A^(L, v) = A^(L, -v 


f(L)- 


Li(3-a) +il3“al j 


( 21 ) 


In L : a = 3 

’ Mcl Iwaln and Fillius (1975) have pointed out that Fig. 1 
shows f(L) — L^ where which is possible in the absence 

of any losses only if a is negative. Since serious synchrotron 
losses yield exponential dependences on M and L, and since the 
weak synchrotron loss criterion is satisfied oy these particles 
within "the limits of experimental measurements of the diffusion 
coefficient, another^ loss mechanism must be active. ! f we drop 
the synchrotron term in equation (1) and replace i t by a pheno- 
menological loss term, then 


^ « |2 _5l n ^ 
bt ^ bL ,2 ^LL bL 


t(L,M) 


f = Q 


( 22 ) 


where t (L, M) is the particle lifetime. 

The inclusion of energy dependences is straightforward 
here; letting and t(L) = where M is 
the first invariant arbitrarily normalized. 
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M+(a-2)L5^.i:J--g ~L f(L; M) = 0 (23) 

bL^ 0^0 

Using the boundary conditions f^L « Lj^ = 0 and f^L = \.^ = 
the solution of (23) is 


and 


f(L, M) = ( 


L>,i(3-o) 


I (xr) I (xb) 


I (xa) I (xb) 

V "V 




T (xa ) n (xb ) 

“V" 

(24) 


a + p -2 


la + p - 21 




S3— ? (m+6 ) 


o o 


-r « L-i(a+3-2) b » L:Ha+0-2) ^ L-'^(a+g-2) ^ 25 ) 


For weak losses, the solution (24) behaves like equation (21). 
For strong losses, the solution behaves asymptotically as 


^ „ ..X . , L , % ( 3-a ) ^ /a sinhrx (b - r) 1 
f (L, M) -V (^) fj^ s±m^^ --Kr\ 


(26) 


Arbitrary values of B yield exponential dependences on L, But, 
for the special case 0 = 2 - a it can be shown that (24) reduces 
to 


T ^(3-et) T %(3’0t)C 
o o 


Q. - (L^/L) 

•[l-(VVi3-a)c^ 


(27) 


where 


4 


(28) 


C = 


1 + 


_^-3s(m+6) 




then outside of the surface absorption region |jL^/L) 


the solution scales like 


f(L, M) « 

o 


(29) 


where 


g = i(3-a) +i|3~a|C 


(30) 


Solving for the magnitude of the particle lifetime 

V,- S 5 ■ ■ 

which is, of course, what will result if we assume a power law 

for f a priori and insert it into (22) . If m + 6 =* 0 and if 

we choose a = 4, 3, 2, then estimating the slope from Fig. 1 

as q- 4, we have as an estimate for D_t = 

O O XO xz 

The fact that the observed distribution behaves like a 
power law in the loss region, L«10, supports the argument that 
t(L) 'vl^~ , This situation, as stated above, is only the require- 
ment that anomalous losses just balance the diffusion source 
and that f behave as a power law in L. . 

If the GorotatiGn drift dominates over the curvature-gradient 
drift at Jupiter for these particles, we ejcpect m- 0 (Goroniti, 
1974) . Figure 1 demonstrates a weak dependence on M inside ' 
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the loss region g7^g(M)> which supports the 

tion 5 =s 0. 

Although the loss mechanism has not been confirmed, loss 
cone precipitation is the leading candidate. For a dipole field 


the minimum precipitation. lifetime on strong pitch-angle dif- 

. 2L^V 

fusion is (Coroniti, 1974) = — -= — r- which is independent 

c(Y^- 1)““ 

of energy for relativistic electrons. For the size of the dif- 
fusion coefficient suggested in section II lO^^^sec”?") , 

we find that for D ® t >> t and that if pitch-angle 
diffusion is the loss mechanism, the weak pitch-angle diffusion 
regime prevails. 

2*ct 

We note that if the particle lifetime presists as r = 
into the synchrotron emission region and that no other, dif- 
ferent mechanism takes over at small L< 2 values, then if m = 5 0, 


synchrotron losses can be included and all of the formulas of 


section II still hold with the only substitution 

, ^ 

1 1 + 


3 -^g , 

+ 11/2 ^ 


tion ..(19) .kSe, of course, still valid. 


The conclusions reached from equa- 


b) Volume Emissivity 

On the assumption that weak synchrotron losses persist 
up to the surface of Jupiter (i.e., > 40 10 ^^sec ^) and 

that the lifetime x = x^L^“^ still applies for L < 3 , we ask 
how does the equatorial emissivity vary with L. If the decay 
in the emissivity for L<l.fa is due not to 1) synchrotron 
degradation or 2) a different, very intense loss mechanism 
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which commences just inside of L<2, then the decay is due to 
surface absorption, if m + 6 *= 0 and c = 


tion (27) becomes 


1 + 


Voj 


f equa- 


f(L. p) = 


4(3-a)4l3-alc4(N+2), 


[l - (1/L)*^'^*^]pN+2 


P 


( 32 ) 


The monochromatic emissivity is defined as 


f* ^ -ii* ^ ^ 

r\fXn) = Jdp Pf(p, n) f(x, p) 


( 33 ) 


where P^(p, n) is the power radiated at the frequency f in the 
direction n by an electron with momentum p. For simplicity 
CCS urns the e I actrons radiate at only the characteristic frequency 


Pf(p, n)®Y^B^(L) 


( 34 ) 


For relativistic particles where equation (3? ) scales 


as 


^f 


.J_r X, 6(f - u) 

(n)-g(L) N-1 f i ( n - 1 ) (35) 

L P J 


and we have 


,.,^(;;)_,_i(3-a)+4l3-alC-3(N+3/2)ji , (,/L)'3-atC^ f-i(K-l) 

( 36 )' 
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The emissivity maximizes due to surface abosrption (no 
synchrotron losses) at L = L' where 


1 


' = 1 + 

3(N + 

u 


(3 - ate 


3/2).-l{3-a)--^l3-alC 


Z 


1 3-a I c 


Til 2q 4- g ~ 3 _~1 
^ 2(nV 3/2)"^ qj 


2q-Ki-3 


( 37 ) 


IV DISCUSSION 


The electrons measured by Pioneers 10 and II were in the 
weak synchrotron loss regime and equation (32) is relevant for 
them. The uni ^di recti onal flux j (x, p) = p f(x, p) and the 
onni-directional integral flux is- then approximately 


4(3-a)4l3-alG4(N+2) 


(38) 


For L>8, where losses are not apparent at high energies, if 

aa4, J(y > r) « ^ . In this region Ni*«3.5 for* electrons 

-8 2 

of energies E > 35 MeV and J(y >70)=-L * , consistent with observati 

If we suppose that the diffusion coefficient is large enough 
to transport the significant synchrotron emitting electrons 
to the surface, and also that the characteristics of the anoma- 
loxis loss mechanism do not suddenly change for L< 2, then the 
electron flux should drop because of absorption at the surface 
and the emissivity should peak at L"* given by (37) . For reason- 
able values of parameters L' is insensitive to q and a: if 
we choose a = 4 and N = 1, L' = 1.13-1.17 for q = 0.5. If 
the resolution of the volume emissivity is high enough to pre- 
clude these values, the particles must undergo serious synchro- 
tron losses before reaching the surface. This conclusion puts 
an upper bound on the magnitude of the diffusion coefficient. 

With the assumption that the 3000 MHz emitting particles dif- 
fuse to L = 1.8 without serious synchrotron degradation, we 
obtained an estimate for the minimum diffusion coefficient of 
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0.88, 2.0, 4.6, 11 < Dq min ^ a = 4, 3, 2, 1. Since 

the 10.4 cm emissivity decreases significantly by L ~ 1.4, we 
can estimate an upper limit on by assuming that the synchro- 
tron loss criterion (20) is not satisfied at L = 1.4. Substi- 
tuting f = 3000 MHz and L = 1.4 in (20) , the range of D values 

w O *N. 

is bounded by - - • . _ ^ 

0.88, 2.0, 4.6. n <0 qX 10^°<4.5, 7.9. 14, 26 (sec“^) 

(30) 

for a = 4, 3, 2, 1. The criteria of minimum synchrotron degrada- 
tion at L s 1,8 and serious losses at L = 1.4 limits the range 
of 0^ to less than a factor of 4, 

Mogro-Campero has reported (see Fllllus et al«, 1975) a 

. > . d 

measurement of the diffusion coefficient at lo of D(L) = D, (■?— r] 

io\5.9/ 

where 1 ,5 x 10'® < Oj^ < 3, 1 x 10”^ sec" ^ and 3.6 < <: 4.0. Using 
the extreme values of gives 0,12, 0. 25 < D^x 10^^ < 2.6, 5.2 
for n^ <= 4, 3.6, We note that in (30) the actual value of 
is likely to be closer to 0^ ntin» since the upper value of M 
(352 MeV/G) which we have used corresponds to the 3000 MHz (— 10.4 cm 
emitting particles, and the 10.4 cm emissivity peaks near L = 1,8. 
Thus, if we estimate 0^ as a probable value for D^, this 
value agrees well with Mog ro- Camper os ' measurement. Simpson 
et a1. (1974) have reported a measurement of the diffusion coef- 
ficient 2x10 ^ < 0|^ < 10'® sec'^ ; assuming a. = 4, 1 .5 < D^x 10^® < 7,7, 
which Is also consistent with equation (39 ). 

Birmingham et a1, (1974), using a pre-Pioneer 10 value 
of 10 Gauss for the magnetic field, concluded that a diffusion • 


10 


•9 1 95 •*! 

coefficient of 1.7x10 L sec gave the best fit to the 

10.4 cm emissivity. We note that equation (20) gives for == 

—1 

and 0=2 the value = 1.7 ^ lo sec when f = 3000 MHz 

o min c 

and L = 1.8. 

Our conclusions can be ordered in decreasing L value: 

..1) The fact that the phase space density obeys a power 
law reasonably well in the loss region 3 < L < 10 of Fig. 1 sug- 
gests the electron lifetime varies as t «= and that the 

relation (31) holds in this region, j 

2) The inclusion of synchrotron losses which are charac- 
terized by exponential dependences in both the particle energy 
spectrum and particle flux, J(L ) , in the strong synchrotron 
loss regime demonstrates that strong synchrotron losses do not 
occur for 3000 MHz emitting electrons for L 5 I. 8 . This pro- 
vides a lower bound for the magnitude of the diffusion coef- 
ficient given by (39) . 

3) If the characteristics of the los^ mechanism do not 
alter and intensify at L£l.8 then the decay of the 10.4 cm 
volume emissivity is concomitant with strong synchrotron losses 
and the diffusion coefficient is bracketed by equation (39) . 

4) The scaling of the particle lifetime as L “ cind the 

g 

magnitude, which is estimated as € x 10 sec from equations (39) 
and (31) for a = 4, indicate that if pitch-angle diffusion of 
the electrons is the culprit, the weak pitch-angle diffusion 
regime prevails. 

5) In the strong synchrotron loss limit the power law 






dependence in energy of the di^'^tribution should collapse and 
behave exponentially as equation (B.8) . An ■ in situ measure- 
ment of the distribution .( > 352 Mev/G) inside L < 1.8 would test 
the hypothesis that serious synchrotron losses are present. 

We note also that for a = 4 and =» 8, 8 x lO'^^sec”^, equation 

(20) requires the 200 MHz (150 cm) volume emissivity to peak 
at Zj ^ 1 • 5 . 
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Figure tap I ton 


Figure 1. Phase space densities reported by Well wain and FMIius 
(197^) from Pioneer 10 data, evaluated at constant 
first invariant M, 
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APPENDIX A 


We must evaluate 


where 


I 




dt e 

l«D+e 


”cp(t) 


(AJ) 


q>(t) 


M\a) r'’(\b} - r"^(xb) 1. (>.aj 

»« • «t %*» •• 


(A. 2) 


If N Is an integer, the integrand is sing I e- valued and 
analytic except for simple poles of cp(t) on the negative real 
axis and the pole at the origin. Also, if N = 1 we must intro- 
duce an upper limit cutoff for to satisfy equation (11). 

If we close the contour by a semicircle to the left, which gives 
no contribution, since at best ep(t)— e^^ compared to e^ (a 
rigorous treatment can be found in Cars law and Jaeger, 1959), 
then by residue theory 



N-1 


(A) [eS(t)j +r(M)£ %i^esl*cp(t )] (A. 3) 


’t=0 


res 


Note that the ant i- symmetric combination cp(t) can be written 


q)(t ) 


l^/Xr) K (Xb) - l,,(Xb) K ,(\r ) 

r(xa) K'"{xbj - l'’(Xb) K^’(Xa) 
V V V V 


(A. 4) 


If 
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T 

e z 


and 


Y (z) = 
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J (z ) cos vTT * J (z ) 

-It 

sin vTT 
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then 


J (zr) Y (zb) - J (zb) Y^/zr) 
. 1>t‘> = /(za) /'(zb)- /(zb) Y TzTy 


(A.5) 


If the roots of -J (x> Y (x^) - J (x^) Y (x) =0 are denoted 

v' V a V a V J 2 

ay , the poles lie in the t plane at t « “ “ ■^v,n* 

The second term in (A, 3) can be expressed as 


" I (xr) I (Xb) - I (Xb) I (Xr) 

= ^ 

nal ‘ 


I (xa) 1 (Xb) - I (Xb) I, (xa)l 

. V ■* V V V , J 


(A. 6) 


t=-t 


v»n 


Using the property that the Wronskian *-vl ^ 

it is Straightforward to show that 


.t^l^(Xa) l.^(xb)- l^(Xb) ..^(xa)]|^^_ 


Sin mtt 
TT 


J^(y b^ - J^fv n®) 

\r' . M . n. / ■ x »_-V\i-.n / 


v,n 


(A. 7) 


and the second term of (A. 3) is evaluate^ as 


N-1 


-2 


I? ="r(N)(-x^) I -2i[^IrrT 

' n=l Yy,n 


J^ilv.,n?) j 4l.VL.n! ’) rj (y b') Y fv r') - (r 

• j2(y b) - j2(y a')'- 

v\ v»n / v\' v,n / (A 8) 
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The first term in (A. 3) Is due to the pole at the origin and 
is 


I, - = 0) (A.9) 

Jls=0 

Now<p(t:) ° utilizing the series expansion of 

’rr we have ' ^ 

±v 

CD n 

D(r, b, t) = y ^ (iX)^” t" 
n=0 k=0 




in (n- k)l r(v + k + 1 } r( * V + n - k + 1 )J 

(A. 10) 

from which we must evaluate (t = 0) . The general result 

\ 

is cumbersome/ so we shall give only the leading and first order 
terms , Define 
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(A. 11) 


A (L) = — b^a/r IV f f ^b l ^ , ^ 
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Then 
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Special cases are 
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l-rVb^-^rVb^ln t./r][ t - a^/b^] 

-■(r-- a)]^ (A, 


As mentioned previously if N = 1 we should introduce a cutoff 
at high momenta fQ(M) » " ' $/ ’z^(J^q ~ modification alters 

the previous calculations only slightly and the result Is, for 
N » 1 
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APPENDIX B 

a) If the arguments of all the Bessel functions in equation? 
(14) and (A.l) are large, the asymptotic expansion of 
produces 


. 2 
a 



In which we have retained the dominant exponentials of both 
numerator and denominator. Inverting the denominator yields 

q)(t ) ~/a/r 





where D 2 ^(z) is the parabolic cylinder function {see Abramowi tz 
and Stegun (1970) for definitions of all functions in this paper) 




we have for the distribution function 
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If N is a positive Integer (B. 4 ) reduces to 
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For large w» } 


2N-2 — 7 ( 1 - a/r 

r" *ERFC(W//2)-r^ ^ ^ =-—r 

[Xr(l;a/r)J^ 


(B.6) 


which demonstrates the exponential dependence on M and *> 
for strong synchrotron losses. Equation (B.4) is useful when 
Xaa1,.for which values the series in equation (A. 18) converges 
very slowly. 


b) The intermediate case of weak losses at L =» L is obtained 

o 

by treating xa«1, Xb » 1 and Xr»l. Expand l^(Xa) in small 
argument form and keep only the leading term. If \» ^0, 

K^Xa) “•jr(v)(xa)'^ dominates. We have then 

® + -Sir ■ i 28 x^r^ J 


(B./; 


which gives 
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ir2^2 
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where 


U = 5(v + 2)-N 
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Xr 

/2 


• Va 






for large values of w 


lr2_2 

n2^(w)-w^“e'® 




9) 


gild f(L) exhibits the same exponential dependence as equation 

(ll»4) when \^r^>4. 
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